Pulse Areas in Multi-Soliton Propagation 
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The prospect of self-consistent propagation of more than two pulses contemporaneously through 
multi-resonant media raises open questions: whether soliton solutions exist, and whether a useful 
generalization of two-level pulse Area can be found. We answer these questions positively for the 
case of four pulses interacting in combined V and A fashion with an idealized pair of atomic D-lines. 
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The fully coherent propagation of light pulses in ab- 
sorbers with a dominant resonant transition was given a 
new foundation, superseding Lorentzian theory, by the 
work of McCall and Hahn n\ in the late 1960s. Their 
theory predicted stable solitonic pulses and introduced a 
new concept, the dimensionless quantum pulse "Area" , 
6{z, t), as the time integral of the pulse envelope £{z, t): 
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is the well-known Rabi frequency and d is the dipolc tran- 
sition matrix element for the resonant transition. The h 
in its denominator indicates that the Rabi frequency has 
no classical h ^ correspondence limit. The McCall- 
Hahn "Area Theorem" establishes that 9 — 2rm iden- 
tifies values for Area that are stable in the propagation 
of pulses whose duration is much shorter than relevant 
homogeneous relaxation times such as the spontaneous 
lifetime [2|. 

The propagation of several fields simultaneously 
through multi-level media with several resonant transi- 
tions opens two questions: whether stable soliton solu- 
tions still exist for any such cases and whether pulse Area 
remains a key parameter. Here we examine these ques- 
tions by extending recent work [^, |3| on three- level media 
to more complex four-level media that can be viewed as 
having either double A or double V form for their levels 
and resonances, as shown in Fig. [TJ 

In the context of Electromagnetically Induced Trans- 
parency (EIT) Q pulse propagation in three-level me- 
dia has been studied extensively [ g , III , [sl, p, and for 
multi- level generalizations of EIT [lO|, [ll|, [ij, [l3| . Trans- 
parency in these systems is usually achieved by hav- 
ing a strong quasi-cw field produce a steady-state dress- 
ing of the medium on one resonance, producing a very 



narrow spectral window for essentially lossless propa- 
gation of a much weaker pulse on the other transi- 
tion. In the opposite regime, where the pulses are wide- 
band and short, the evolution of both pulses is coher- 
ently dynamic, leading to soliton formation. This situa- 
tion has been examined for particular multi-level media 

i, i, Q 0, El, [13, E, HlH, El, El, m, H, but the 

role of pulse Area has been largely overlooked. 

Under the condition of pulse-matching [1, , simul- 
ton pulses |l4ll are known to obey a two-pulse Area The- 
orem [l^, l2l| in y-type media and multi-level generaliza- 
tions exist for certain systems 17[. The recent result of 
Clader and Eberly [3| gave an exact analytic solution for 
A-type media in which the pulses are temporally matched 
but have different spatial propagation. Eberly and Ko- 
zlov have shown that the propagation of unmatched two- 
photon-resonant pulse pairs in A-type media depends on 
the so-called dark Rabi frequency [20| and obeys a dark 
Area Theorem l22jl. 
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FIG. 1: A four-level atom interacting with four separate lasers 
with equal detunings. A, and four laser pulses indicated by 
their Rabi frequencies. In the short-pulse regime, we can ne- 
glect relaxation processes, as indicated by the crosses in the 
rates 73 and 74. The sketch can be interpreted as showing 
idealized alkali D line transitions, where levels 3 and 4 are 
P\/2 and P'i/2 manifolds, and levels 1 and 2 correspond to 
F=l and F=2 hyperfine levels of 51/2. 



2 



In this work, we examine four-pulse propagation under 
conditions suggested by Fig. [Tl such that the durations 
of all pulses are short enough to neglect homogeneous re- 
laxation. Hioe has shown 17] that this level scheme per- 
mits simulton solutions (multiple pulses travelling with 
the same group velocity). Here we employ the Backlund 
Transformation method of Park and Shin It 
is robust enough to accommodate both inhomogeneous 
broadening, which we take into account, and more than 
a single pulse velocity. As is evident in Fig. [TJ atoms 
in the four-level configuration can be viewed as a dou- 
ble A- or double t^-type medium and our solution has 
attributes of known A and V solutions. We have discov- 
ered that pulse Area remains a parameter important to 
the evolution of our system and that three "total" pulse 
Areas are constant. 

For the four short optical pulses, we assume that 
the laser fields are linearly polarized and can be 
written in envelope-carrier wave form as Ea{x,t) — 
fa(x, t)e(''^°^~"''*) -I- c.c, where Ea{x,t) is the slowly- 
varying field envelope function, ka is the wave number 
and uja is the field frequency addressing the 1-3 transition 
labelled Similar notation applies to the other laser 
fields and each is assumed to address a single atomic tran- 
sition. Under the Rotating Wave Approximation (RWA) , 
the Hamiltonian for this atom-laser system can be writ- 
ten (in the |1) , |2) , |3) , |4) basis indicated in Fig.[T]) as 
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(3) 



where the Rabi frequency fia governs transitions between 
levels 1 and 3 through the atomic dipolc moment operator 
da associated with the a transition, and the other Rabi 
frequencies are similarly defined. 

Transitions between levels 1 and 2, as well as between 
levels 3 and 4, are assumed to be dipole forbidden due 
to the parity of the levels. We require the detuning to 
be the same for all transitions: A = ^3 — {u}\ + oja) = 

W3 — -I- Wc) = W4 — (wi -I- ijJh) — L04: — (1^2 + i^d)- As 

indicated in Fig. [1] we ignore homogeneous relaxation 
rates, such as 73 and 74. The atomic density matrix 
then satisfies the von Neumann equation 



ih 



'dt 



[H,p]. 
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The propagation of each field is governed by Maxwell's 
equations which, under the Slowly- Varying Envelope Ap- 
proximation (SVEA), become: 



dZ 
'dZ 
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'i^J-b (pii) , (5a) 
-i^d (P24) , (5b) 



where we use the retarded coordinates T = t ~ x/c and 
Z = x/c. The atom-field coupling parameter is fia = 
Mdl^uja/heo, where M is the density of atoms. Similar 
notation applies to /if,, /ic, and pa- The brackets denote 
an average over all atoms to account for inhomogeneous 
broadening. We assume an inhomogeneous lifetime T2* 
and line-center tuning. The distribution function is then 
F(A) = (T2*/V2^)ea;p[-(AT2*)V2], and the averaging is 
performed by integrating over the detunings. The form 
of the von-Neumann equation remains the same in the 
retarded coordinates but the derivative must be taken 
with respect to the retarded time. 

We apply the Backlund Transformation method of 
[3, [IBl to our four-level system by first recasting Eqn. 
^ in terms of the Hamiltonian. This is easily done by 
introducing a constant matrix 



W ■ 



/O 0\ 
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so that the slowly- varying Maxwell equations can be writ- 
ten in commutator form: 



dH 
'dZ 



(7) 



Note that we have required the atom-field coupling for 
each transition to be the same: p = pa = Pb = Pc = Pd- 

Written in this way, Eqns. ([5]) and ([7]) differ from the 
A-system case only in the definition of W. The solu- 
tion method outlined by Clader and Eberly, based on the 
Backlund Transformation of Park and Shin 20] , does not 
employ the explicit form of W and can thus be used to 
generate solutions to our system. 

Given any "seed" soliton solution to a set of non-linear 
partial differential equations, the Backlund Transforma- 
tion generates a higher-order soliton solution from it. An 
obvious exact solution to Eqns. ([5]) and ([7|) is given by 
zero applied fields: 

r!,(r, z) - n,{T, z) = r!,(r, z) = n^ir, z) = o, 

and constant level populations, as specified by this mixed 
state density matrix: 



p(r,z) = p(o,o) 



(\af 0\ 

|/3|2 



\ 0/ 



(8) 



where lap 



= 1. To be specific, we will assume in 



the present discussion that \a\ 



> 



We have applied 



the Backlund Transformation method to the seed solu- 
tion above, which leads to exact solutions for the Rabi 
frequencies and atomic density matrix elements. We will 
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focus on the Rabi frequencies and will not list the density 
matrix elements here. The Rabi frequencies are: 

naiZ,T) = sinusin(/)(Z)f7(Z,T), (9a) 

flb{Z,T) = cosusin^(Z)f](Z,r), (9b) 

nc{Z,T) = smucos(j){Z)n{Z,T), (9c) 

VLd{Z,T) = cosucoB(j){Z)n{Z,T), (9d) 

where the fundamental Rabi frequency f2 has the com- 
plicated form: 



n{Z,T) 



(4/t)^1 + exp [2(^2 - /32)^-^ 



2 cosh [a^nZ - + exp [( 



2/32)kZ- 



and tan(/)(Z) = exp[(/32 ~ q;2)kZ]. 

In the expressions above, r is the nominal pulse width 
and the length scale is set by k, defined as the inhomo- 
geneous average n — (^/2t)(1/(A2 + l/r^)). The pa- 
rameter u is a constant of integration that specifies the 
initial conditions of the pulses. For A-type media, this 
constant is of little physical interest because it merely 
shifts the origin in the spatially infinite medium. In the 
four-level system, however, u has a natural interpretation 
and allows freedom in the choice of the relative pulse am- 
plitudes. 

The pulse propagation behavior is difhcult to infer 
from the complicated formulas above. However, plots 
show that the evolution can be broken down into three 
distinct regimes: regime I, where Z and T are large and 
negative; regime II, where Z,T ^ 0; regime III where Z 
and T are large and positive. 

In the asymptotic limit of regime I, where 
~kZ, —T/t 3> 1, pulses ^la and have sech form, with 
Areas equal to 27rsinit and 27rcosu, respectively, and 
pulses fic and fi^ approach zero amplitude. In regime 
II the complicated 4-level interaction process begins to 
convert pulses and Vtb into pulses f2c and fid, which 
eventually become sech pulses in regime III, where 
kZ,T/t 3> 1. They asymptotically approach Areas 
equal to 27rsinM and 27rcosu, respectively. The input 
and output pulse pairs are not identical, however, if 
a ^ (3. Then the input group velocity Vg ~ c{l + a'^KT)^^ 
shared by U,a and 17;, differs from the output group 
velocity Vg — c(l -I- fP k,t)~^ shared by Sic and f2d, and 
the output velocity is greater under our assumption 
that a > (3. The conversion from pulses Q.a and Q.b to 
pulses Q.C and Q.d is a feature not present in the pure 
"simulton" pulses of Hioe [l3| , which all travel with the 
same group velocity because a = /? is imposed in that 
case. 

We need to probe the stability of these solutions. That 
is, we need to know if these analytic solutions to ideal- 
ized medium and pulse conditions can provide useful ex- 
perimental insights. Thus, we have employed numerical 
methods to calculate the expected behavior of non-ideal 
input pulses in finite media. 



A particular example is shown in Fig. [2] The simula- 
tion was performed with square-Gaussian input pulses 
(i.e., having quartic rather than quadratic exponents) 
with total Area 6t less than the predicted stable soli- 
tonic value of 27r. In contrast to the analytic conditions 
of an infinite medium, i.e., zero amplitude for pulses f2c 
and rjrf in the asymptotic input regime, we alloted a small 
initial amplitude to pulses f2c and Sl^ before they enter 
the medium. 
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FIG. 2: Snapshots of time-matched pulses represented by 
their Rabi frequencies generated from numerical simulations. 
The position of the medium is represented by the vertical 
lines. All pulses are initially of square- Gaussian shape with 
initial areas approximately those of regime I . We have chosen 
Oa = lAn, Ob = O.Qtt, Oc = 0.0027r, and 6^ = O.OOlvr so that 
6t is less than 27r in Frame 1. 



Upon entry into the medium, pulses Via and U.^ are 
quickly reshaped toward the preferred solitonic sech form 
(Fig. [21 Frame 2). The sech pulses travel at a reduced 
group velocity in the medium, which results in their re- 
duced width along the coordinate axis (Fig. [2l Frame 
3). With further penetration into the medium, the weak 
pulses r^c and U,d are amplified in the pulse transfer 
regime as ^la and fib are depleted (Fig. [51 Frames 4,5). 
Finally (Fig. [21 Frame 6), fie and are of stable sech 
form and propagate without attenuation through the re- 
mainder of the medium. The greater width of the "out- 
put" pulses and fi^, in comparison with f2a and fib, 
along the coordinate axis is a consequence of their greater 
group velocity in the medium. 

The parameter u determines the maximum amplitudes 
of the pulses and clearly shows that the maximum am- 
plitude of f2c is determined by the maximum amplitude 
of ^la and, similarly, that the growth of f2d is limited by 
VLh- This behavior can be best understood when viewed 
in the context of pulse Area, as follows. 

Pulse transfer processes were previously observed for 
A-type media by Clader and Eberly [3, |^ . They found 
that the individual pulse Areas, defined by 9{Z) = 
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FIG. 3: Pulse areas in units of tt as a function of propagation 
distance in the medium, indicated by the vertical lines, from 
numerical simulation. The dotted lines show the "total" pulse 
areas 6i, 62, and 9t, as indicated on the right side of the figure. 



shown that pulse Area remains a powerful nonlinear con- 
straint on the evolution of pulses in four-level media, and 
more general combinations of the pulse Areas are stable: 
three "total" pulse Areas remaining constant throughout 
the medium. Numerical simulations have been used to 
confirm the utility of our analytic solutions under non- 
ideal conditions for which such solutions arc not avail- 
able. 
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J^oo ^(-^i T)dT , were constant in the asymptotic regimes 
I and III and that, surprisingly, a "total" pulse Area was 
constant throughout the entire medium. We have found 
that similar behavior occurs in a four-level medium. We 
can calculate the Area of each pulse from Eqn. ^ and 
find 

OaiZ) = 27rsinusin(?!)(Z), (10a) 
eb{Z) = 27rcos'usin0(Z), (10b) 
edZ) = 27rsinucos(^(Z), (10c) 
OdiZ) = 271 cos u cos (j>{Z). (lOd) 

We note that particular combinations of these Areas 
are constant throughout propagation. In particular, we 
find 

0i{Z) = VKiZW + MZW - 27rsinu, (11a) 

02iZ)= y^iZj^TMZW --27TCOSU, (lib) 

0T{z) = ^Mz)\^ + \e2{z)\^^2TT. (11c) 

As in the A-system, the "total" pulse Area 9t is ex- 
actly 27r throughout propagation, for any value of Z. In 
addition, the four-level system has two more relations, 
Eqns. (|llap and (jllbp . restricting the pulse Areas. These 
clearly show that amplification of flc is a consequence of 
the depletion of ila and that, similarly, the Area of il^ 
grows at the expense of the Area of r2{,. 

To check these predictions derived for ideal pulses and 
media we show the pulse Areas corresponding to the pulse 
evolution in Fig. ^ also computed numerically. They are 
plotted in Fig. [31 One easily sees that the total pulse 
Areas grow to the analytically predicted values. 

In this note we have presented new exact analytic solu- 
tions to a four-pulse propagation problem, extending pre- 
vious simulton results into a regime that can be identified 
as four- wave mixing of solitons. The solutions have been 
shown to combine elements of known solutions for V- and 
A-type media. The three regime pulse-transfer behavior 
of A systems identified in [3[ is again present. We have 
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